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Abstract We present early results of the application of 
a method which uses multicolor photometry and spec- 
troscopy for £ discrimination. This method has been 
successfully applied to the pulsating hot subdwarf Bal- 
loon 090100001. Here wc apply the method to QQVir 
(PG1325+101). This star was observed spectroscop- 
ically and photometrically in 2008. Details on spec- 
troscopy can be found in Telting et al. (this volume) 
while photometry and preliminary results on £ discrim- 
ination are provided here. The main aim of this work 
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was to compare the value of the £ parameter derived 
for the main mode in QQVir to previously published 
values derived by using different methods. 

Keywords hot subdwarfs; pulsating stars 



1 The method 

The method we used for mode degree £ discrimination 
has been already described and the results of its ap- 
plication to the main sequence stars presented in sev- 
eral papers. This method is very challenging, since it 
requires observations performed in at least three fil- 
ters. Moreover, time-series spectroscopy could also be 
employed. For this reason, the best objects to apply 
this method to are bright stars with relatively long pe- 
riods and relatively high amplitudes, so they are ac- 
cessible even with small size telescopes. It is possi- 
ble then to achieve good sampling of the brightness 
changes to derive amplitudes and periods of the de- 
tected modes. Although hot subdwarfs (sdB) are rel- 
atively faint objects, in 2004 independent photometric 
and spectroscopic observations of the brightest pulsat- 
ing sdB star, Balloon 090910001 (B al09), allowed a suc - 



cessful application of the method ([Baran et al] 12008) 
For the sam e star, but with high precision multicolor 



photometry, ICharpinet et al] (|2008|) also derived good 



discrimination for a few modes with highest amplitudes. 
From these two papers we can draw one important con- 
clusion. The method may give positive discrimination 
of the £ parameter either having high precision three - 
filter photometry, or a few (at least two) of low precision 
photometry supported by spectroscopy. Encouraged by 
the successful application of the method for Bal09 wc 
decided to try with another pulsating subdwarf B star, 
QQVir. 
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The main goal of performing another observation of 
QQVir was to decide which (or if any) value of the 
I parameter for the dominant m ode, obtained in pre- 



vious works, is the correct one. iTelting fc 0stensen 



( 20041 ) analyzed spectroscopic data and stated that 
the main mode is consistent with £ = 0. Quite inter- 



estingly, ICharpinet et al.l (|2006), based on photometry 
data only, derived I = 2 from frequency matching. The 
Charpinet et al. ( 20061) result is even more surprising 
as it indicates unusually high intrinsic amplitude of the 
mode with 1=2 mode. It is interesting then to use 
both kinds of data, photometry and spectroscopy, to 
constrain the value of the mode degree for the domi- 
nant mode in QQVir and if possible, either to confirm 
one of the already obtained values or to exclude both 
of them. 



2 Observations 

From the previous observations we knew that the pul- 
sational period of the dominant mode is relatively short 
(~2min.). Along with a low brightness (~13.5mag in 
the V filter), data required for the method were not 
easy to obtain, particularly with small size telescopes. 
For this reason, we performed only one filter photom- 
etry on each telescope with the exposure time short 
enough to get appropriate sampling but not necessar- 
ily long enough to have high precision data. Details 
on the observations are presented in TablcHJ while all 
photometry gathered in this project is plotted in FigJTJ 
This figure also includes spectroscopic data which have 
been used in this method. Details on sp e ctrosc opy can 



(120101 ) in this 



be found in the paper by ITelting et al. 
volume. 

Data in different filters were taken at different tele- 
scopes. This is the reason why there are not the same 
number of nights for all filters. With respect to the 
spectrum of the star and spectral response of the CCD 
we used, the highest signal has been obtained in the 
V filter. These data have been used for precise fre- 
quency determination of the dominant mode. That was 
also the reason why observations in the V filter have 
been performed for the longest time from the sites with 
easy access, like: Suhora, Mercator and Baker. In ad- 
dition to the V filter, data in the U filter were taken 
at Loiano and Nordic Optical Telescope (NOT) during 
three nights in total. Fortunately, the dominant mode 
is clearly resolved with only one night of data. Some- 
time later, it turned out that, quite luckily, almost at 
the same time as our campaign, absolutely independent 
data in B filter were taken during the course of a dif- 
ferent project. That enables us to use three filter data. 
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Fig. 1 Light curve of QQVir in three filters. Data in 
different filters are colored according to the bandpass and 
denoted on the left side. Red data are for the radial velocity 
curve and the scale on Y-axis is not relevant for these data 
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Fig. 2 Amplitude spectrum of QQVir in V filter. The 
dominant mode is clearly seen at ~7.25 mHz. Significant 
signal is also present at higher frequencies. This spectrum is 
constrained to the region where the dominant signal appears 
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The B filter data were not crucial, but this method of 
mode discrimination improves with the number of fil- 
ters used. 

Although the observations in different filters were 
not conducted simultaneously, we can use published 
results to see if the amplitude of the dominant fre- 
quency changes wi t h tim e. Results obtained b 



Table 1 Observational log for QQ Vir data. All data were 
taken in 2008. 



Telting fc 0stensen I ( 2004 ) and ISilvotti et all ( 200 



0) 



have radial velocity and B filter amplitudes consis- 
tent with ours. Those observations were performed 
independently and at differing times during 2003. We 
therefore conclude that the amplitudes are not chang- 
ing significantly and remain stable to within the er- 
ro rs. Moreover, result s from the multisite campaign 
of ISilvotti et al.l (|2006t ) indicate there are no nearby 
frequencies to the dominant one, even at very low am- 
plitudes. This suggests that our data do not suffer from 
amplitude variability or the effects of nearby frequen- 
cies. 



3 Fourier analysis 

A prewhitening procedure was applied to derive fre- 
quencies, amplitudes and phases of the pulsation 
modes. First, we assumed that amplitudes can vary 
linearly. As it turned out, the change was insignificant 
(within the errors) and we have decided to assume con- 
stant amplitude and redid the prewhitening process. 
Phases were calculated relative to epoch 2454552.0 
which was chosen arbitrarily. Amplitudes and phases 
for the dominant mode at frequency 7.2554776 ± 
0.0000004 mHz (derived from V filter data) derived 
from the data in U B V filters, as well as from the radial 
velocity curve are presented in Tabled A raw ampli- 
tude spectrum calculated from the data in V filter is 
shown in Figj2] It has been limited to the region where 
the signal appeared, which is 6-9mHz. The amplitude 
spectra calculated from other filters looks similar, how- 
ever the noise level and resolution is worse as less data 
were obtained. Since the main aim of this work is to 
derive the most likely value of the dominant mode, the 
other modes are not included in Table[3J Nevertheless, 
during the prewhitening process they were taken into 
account to check if their removal has any influence on 
the frequency and/or amplitude of the mode in ques- 
tion. If so, they were included in the solution. 



4 Application of the method 

Having derived frequencies, amplitudes and phases we 
could apply the method to discriminate i parameters 
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for the dominant modes detected both in the radial- 
velocity and brightness variations. Discrimination of 
£, does not depend on the observations alone. It also 
depends on the model parameters, in particular effec- 
tive temperature, surface gravity, flux and limb dark- 
ening derivatives over effective temperature and surface 
gravity. For QQ Vir we assumed t he effective tempera- 



ture a nd surface gravity derived bv lTelting k, 0stensen 
( 20041 ). T cff = 34800 K, log 5 = 5.81 dex. Then we 



needed to calculate flux distributions in the range be- 
tween 300 and 700 nm for a grid of models with T c ff 
ranging from 32 500 to 37500 K and logg, from 5.5 to 
6.0 dex. To do this, we used a grid of models prepared 
for the interested range in T e g- and logg. The resulting 
flux distributions were multiplied by transmission func- 



■1 



tions of the filters to get integrated fluxes in the three 
photometric bands - U BV - that we used. 

For the same grid of models, we also calculated the 
specific intensities as a function of fi = cos 9, where 9 
is the angle between the line of sight and normal to the 
stellar surface in a given point of the stellar disc. Using 
these intensities, the limb darkening law was derived 
fitting the coefficients c and d of a function 

/(//)= /(0)[l-c(l- M )-d(l-VM)] (i) 

which was found to reproduce the calculated changes 
of I(fi) best. This was done for different combinations 
of input parameters of the model, thus allowing the 
calculation of the necessary derivatives. 



5 Discussion and conclusions 

A disadvantage of the method we applied for mode de- 
gree discrimination is that we cannot derive a value of 
£ in a direct way. We must use an iterative process and 
employ a merit function to discriminate which mode is 
the most probable. The minimum of the function (if 
it exists) indicates the possible value. In this study we 
applied this method to discriminate between divergent 
solutions for the main frequency of QQ Vir. Results of 
our calculation are shown in Fig(31 The results clearly 
show that the best fit was obtained for £ >0. However, 
the minimum of the merit function is not very deep and 
cannot exclude i=l and 1=2. It can, however,, exclude 
I >2. This does not help to achieve the main goal of 
unique mode identification, but is consistent wi t h pre 
vious work. 



In the case of Bal09, iBaran et al 
X 



excluded any £ value if the y 2 function was at least 3 



times higher than the smallest value of this function. If 
we follow this rule we can also exclude the I > 2 value. 
On the other hand, we can use the probability Q that 
the merit function for the correct model will exceed a 
given value by chance. The Q values for the £=0, 1 and 
2 are 0.7, 0.4, 0.13, respectively In this case none of 
the first three Ps can be excluded. However the chance 
that £ = is wrong is more than 5 times smaller than 
for 1 = 2. 
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Fig. 3 Results of application of the method for mode de- 
gree discrimination. The plot shows \ 2 as a function of the 
I parameter. In the bottom of the figure the precise values 
of x 2 for each I are given 
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Table 2 Results of the prewhitening process. Only results for the main mode at 7.2554776 mHz are presented. Phases are 
given for epoch 2454552.0 
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